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RELATIVISTIC ELECTROMAGNETIC WAVES IN AN ELECTRON-ION PLASMA 


Abraham C.-L. Chian* and Charles F. Kennel 

Institute of Geophysics and Planetary Physics 

University of California 

Los Angeles, California 90024, USA 

ABSTRACT 


High power laser beams can drive plasma particles to relativ- 
istic energies. An accurate description of strong waves requires 
the inclusion of ion dynamics in the analysis. The equations govern- 
ing the propagation cf relativistic electromagnetic waves in a cold 
electron-ion plasma can be reduced to two equations expressing con- 
servation of energy-momentum of the system. The two conservation 
constants are functions of the plasma stream velocity, the wave 
velocity, the wave amplitude and the electron-ion mass ratio. The^~ 
dynamic parameter, expressing electron- ion momentum conservation in 
the laboratory frame, can. be regarded as an adjustable quantity, a 
suitable choice of which will yield self-consistent solutions when 
other plasma parameters have been specified. Circularly polarized, 
electromagnetic waves and electrostatic plasma waves are used as 
illustrations*. 


INTRODUCTION 

Recently, there has been considerable interest in the subject of 
strong electromagnetic waves in plasma because of its 'amplications; 
to laser-plasma .interaction and pulsar electrodynamics^ ' • 
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In the presence of an intense laser field, the plasma particles t 
can acquire quiver oscillation energies exceeding particle rest mass 
energies such that relativistic effects of particle mass variation 
have to be taken into account^. Several studies have demonstrated 
that relativistic effects are of significant importance in such 

rption^^ , soliton forma- 
ile modification^ ' and pro- 

In earlier works on strong traveling waves (see e.g. Kef ♦ 10-12) 
only the electron motion was considered, whereas the motion of posi- 
tive ions was ignored. The results obtained thus apply to compara- 
tively small laser intensities for which the ion dynamics is negli- 
gible because of its large rest mass. Nevertheless, when a very in- 
tense laser field is present, ions can attain considerable quiver 
velocities. For instance, it has been suggested that the ion dyna- 
mics may be important in the later stages of the resonant absorption 
process^ , in a recent study of the inverse Faraday effect it 

was shown that the DC magnetic field induced by a circularly polar- 
ized laser can be greatly reduced by the effects of ion motion. 

Therefore, an accurate description of strong waves requires the in- 
clusion of ion dynamics. 

In the theory of relativistic electromagnetic .waves in an elec- 
tron- ion plasma a dynamic parameter, relating electron and ion mo- 
menta, necessarily appears. In previous papers (see e.g. Ref. 14)., 
strong wave solutions were obtained by treating the dynamic para- 
meter either as an arbitrary constant or as a plasma drift velocity. 

This leads to mathematically correct, but physically misleading 
solutions. 

The purpose of this paper is to elucidate the role played by 
the dynamic parameter in the determination of self-consistent strong 
wave solutions. It is shown that, the problem can be reduced to two 
conservation relations for the energy-momentum of the system. In 
the laboratory frame, the dynamic parameter expresses conservation of 
electron-ion momenta. In addition, there is an amplitude parameter 
that expresses conservation of electromagnetic field-particle energy 
density. These two parameters are functions of the plasma stream 
velocity, the wave velocity, the wave amplitude and the electron- 
ion mass ratio.. We will conclude that the dynamic parameter must be 
regarded as an .adjustable quantity, equivalent to an eigenvalue of 
the system, a suitable choice of u^hich will yield self —consistent 
wave solutions when the value of the amplitude parameter and other 
plasma parameters have been allocated. 


laser-plasma phenomena as resonant abso 
tion^, self-f ocusing^ 7 ' , density prof 
duction of large DC magnetic field ^ ' • ' 
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CONSERVATION RELATIONS IN THE LABORATORY FRAME 

Consider a cold two-fluid electron- ion plasma. The laboratory 
frame is chosen to be an arbitrary frame in which a traveling wave 
propagates with a constant velocity cz/n [where c - speed of light, 
n = index of refraction and 2 = (o,o,l)], the plasma moves with a 
certain stream velocity (where -c < V g < c) , and the solutions 
depend on space and time coordinates only through the combination 
0 = t -nz/c. 


The governing equations written in terms of the phase G are the 
relativistic eauations of motion 
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where y = (1 - v /c ) and a = (e,i). 
a a * 

Conditions describing the average particle number density and 
flux can be obtained by taking the phase-average of Eqs. (3) and 
(6),. which gives 
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where the angular bracket denotes averaging over one period m 6. 
In nonlinear theory, it is convenient to define the plasma stream 
velocity as the ratio of the average particle flux to the average 
number density( lTl 5) . It, follows from Eqs. (7) and (8) that the 
electron and ion streaming velocities are equal 
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A stationary plasma corresponds to' the particular case in which 


It is easy to show that -V s , as defined, by Eq . (9), Lorentz 
transforms like any velocity. Suppose there are two frames S and 
S, where S has a velocity Vz relative to S’. Then, the averaged 
particle flux and number densities in two frames are related by 

<Nv.> = <N'v,’> , <Nv > = T(<N'v '> - V<N’>) ; (10) 

-L — 2 ^ 


<N> = r(<N’> - V<N'v ’>/.c^) > (11) 

where ^ = (x,y) and T = (1 - v 2 /c 2 )~ h . Note that phase averaging is 
a Lorentz- invariant operation. Dividing Eq. (10) by Eq. (li) gives 
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which indeed satisfies the Lorentz transformation for velocities.. 
This suggests that if in S' the plasma has a certain stream velocity 
VI, then an observer in S that £as a velocity V* relative to S will 
"see" a stationary plasma with V = 0. Hence, the nonlinear i-sper- 
sion relation for a streaming plasma can be obtained from the non- 
linear dispersion' relation for a stationary plasma by a Lorentz 
transformation^ ^ • 


The behavior of electron and ion number densities follows upon 
integrating Eq. (2), giving 

* 

N = - N - / - 
a 1 - nv /c 
az 

where N* = <N > (1 - nV /c) . 
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An integration of Eq. (5) gives 


B = nz x E -f <B> 


( 14 ) 


where <B> denotes the magnetostatic field. In this paper an unmag- 
netized plasma is considered so <B> =0. 


A conservation relation for electron and ion momenta can be ob- 
tained from the transverse and longitudinal components of Eq. (1), 


namely ^ ^ 


u. t + pu . = D, 

i- e- - 
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where u = yv/c, p = m / m. and D is a constant vector. Eqs*. (15) and 
(16) can be combined into a single equation 


u. + pu - n(y. + py )z = D 
i e i e 


(17) 


Thus, the dynamic parameter D expresses conservation of electron- 

ion momenta. In the absence of a wave, v = v . = V , Eq . (1 7) be- 

* e i s* M 

comes 

(1 + V) (1 - V 2 /c 2 ) _ls V - n(l + u) (1 - V 2 /c 2 )~ h z = (18) 

s s s 

which shows that D is a function of V s _^ n and p. It will be seen 
later that for large wave amplitudes D also depends on the wave in^ 
tensity. In linear theory, = V g + ^p a ( where ^denotes small 
perturbations), Eq. (17) reduces to 


v . 4* pv =0 

ip ep 


(19) 


Hence, in this limit, D is still given by Eq. (18) because the aver- 
age values of the perturbed particle velocities are zero. 


An energy conservation relation can be obtained as follows. 
A scalar product of Eq. (1) with v a /c^ gives 
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where Eq. (13) has been applied. Adding Eqs. (20) and (21), then 
making use of Eqs. (5) and (6), it yields 
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An integration of Eq. (22) then yields a conservation relation for 
the electromagnetic field-particle energy density 
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where W is a scalar constant. The energy constant W can be consid- 
ered an amplitude parameter that characterizes the magnitude of wave 
intensity, as will be shown later. In the absence of a wave, Eq. (23) 
reduces to 


W = (1 + i/p) (1.- V s 2 /c 2 ) ~ h 


(24) 


W must exceed the above value for a wave solution to exist. 


CIRCULARLY POLARIZED ELECTROMAGNETIC WAVES 


Purely transverse, circularly polarized electromagnetic waveS 
have been studied extensively (see e.g. Ref. 1, 17—19). For the case 
of an unmagnetized electron-ion plasma str eam ing in the wave direc- 
tion, the dispersion relation in the laboratory frame is. 
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where co_ a = 47T<N G >e /m a and v = eEmax^p 1 ^ 0 an i nvar i ant parameter 
that measures the wave amplitude. Eq. (25) relates the wave velocity 
the plasma stream velocity, the electron-ion mass ratio and the 
wave amplitude. The aim of this paper is to demonstrate that similar 
dispersion relations for other, strong waves can also be obtained^ 
through a proper treatment of the conservation parameters WAnd D. 


ELECTROSTATIC PLASMA WAVES 

As an illustration, consider the case of electrostatic plasma 
waves. For longitudinal waves, u a = (0,0, u a ) and the two conser- 
vation relations Eqs. (17) and (23) reduce, respectively, to 
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where t = w pe 0 (1 - nV c /c) . Evidently, Eq. (27) indicates that an 
oscillatory solution exists, with a period of oscillation given by 
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where the turning points u ^ , 2 are determined by the equation 
Y e + = W. Note that is related to u e and y e through Eq. 

(26). The general dispersion relation (28) recovers the electron 
plasma result P h e ^ on term Y^/p is ignored. 

Before discussing the general solutions of Eq . (28), consider 
first the phenomenon of wavebreaking which may occur for subluminous 
waves and is relevant to the laser-plasma interaction^^. A condi- 
tion for the occurrence of wavebreaking in the laboratory frame can 
be obtained as follows. First notice from Eq. (13) that, for a 
stationary plasma, a traveling wave solution exists (i.e. 0<N^<«>) 
only if the condition 


v < c/n 
a 


(29) 


is satisfied. Upon transformation from the stationary plasma case 
to the streaming plasma case^^ the condition (29) then becomes 


where 
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Hence, for subluminous longitudinal waves, the electron and ion 
velocities are subject to the conditions (29) and (30), whose 
violation implies the occurrence of wavebreaking. A graphical dis- 
play of Eq. (30) in Fig.. 1 shows that, for a given wave velocity. 


Normalized Critical Particle Velocity 
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Figure 1 . 


Variation of v c 
for V s /c = 0.5. 


/c and 1/n with normalized wave velocity 
Shaded region indicates wavebreaking. 


the ran-e of v„/c permitted for traveling wave solutions is bounded 
by th, 1/n cur?! and the y c /c curve. The shaded regron represents 
the domain, in which wavebreaking takes place. 


The problem of wavebreaking can further be clarified b> ex am in 
Lng the conservation relations in the wave frame, for which there is 
a/ time dependence. In the wave frame, the governing equati 

tome 


d f \ a r 
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which shows that oscillatory solutions exist only if lies in the 
interval 

W < W < W (38) 

mm max 

J 

where V - =. -[D, 2 - (1 + v) 2 ]^/p and W max = -[(v “ Dj) ~ O'Vu- 

It is easy to see from Fig. 2 that wavebreaking occurs if exceeds 
w max . Furthermore, it shows that Wj determines the amplitude of the 
wave oscillations. 


The general properties of electrostatic plasma waves in the 
laboratory frame can be computed from Eq . (28). When evaluating a 

particular solution it is necessary to decide what values should be 
assigned to the two conservation parameters W and D. Our wave frame 
analysis indicates' that W can be considered as an amplitude para- 
meter. To obtain a specific solution with given plasma stream velo- 
city, wave velocity, the electron— ion mass ratio and the wove ampli- 
tude, the dynamic parameter D must be adjusted accordingly so as to 
render self-consistency to the solution. Fig. 3 shows the variation 
of the cut-off frequency u)/u) (where = w pe 2 + w pi an ^ n -0) 

with v c /c (where v a = |vj - V 2 I, v l>2 being the turn:i - n g points of v e 
oscillations) for different values of plasma stream velocity. The 
computations are carried out for a hydrogen plasma (y = 1/1837). 

In Fig. 4, the variation of the self-consistent dynamic parameter D 
with u a (where u a « |u x - u 2 |) for a stationary plasma is displayed 
for different values of n. It confirms that the dynamic parameter 
is a function of wave amplitude. Note that for small wave ampli- 
tudes, D stays very close- to- the value giv$n by Eq . (18). 


CONSERVATION .RELATIONS IN THE SPACE- INDEPENDENT FRAME 


For super luminous electromagnetic waves the analysis can be 
simplified by referring to the space-independent frame which has a 
velocity ncz with respect to the laboratory frame. The basic equa- 
tions in this frame are 
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The above system of equations can be combined to yield two conser 



for indicated values ’of n; 


Figure 4. Variation of D with u 
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Eq. (43) corresponds to Eq . (17) in the laboratory frame. Upon 
transforming to the laboratory frame, the transverse components of 
Eq. (43) gives Eq. (15) with D 9 ^ = , while the longitudinal compo- 

nent of Eq. (43) gives Eq . (167 with I>2 Z - (1 - n^)“"^D z . Hence the 
dynamic parameter expresses conservation of electron— ion momentum in 
the space-independent frame. Eq. (44) corresponds to Eq . (25) in 
the laboratory frame and can be rewritten as 

du \ 2 

3T +Y. + Y«/P-W, . < 45 > 





where T 2 = w pe2^’ ^peZ = AttNc /m e end W 2 — (1 — n )^W» Therefore* 
the amplitude parameter expresses conservation of electromagnetic 
field-particle energy density in the space-independent frame. x For 
the special case of Bg = D 2 Z, W 2 must exceed [d 2 2 + (1 + .v) 2 ]^/v for 
a solution to existC 1 ^). . An elegant formulation of the energy- 
momentum conservation relations in terms of the Maxwell r s stress 
tensor in the space-independent frame can be found in Ref. 21. 


The two energy-momentum relations (A3) and (A5) can be used to 
study all possible wave solutions. As in the case of electrostatic 
plasma oscillations, a proper handling .of the conservation constants, 
§2 and W 2 , is essential for obtaining self-consistent solutions. 


The case of circularly polarized electromagnetic waves serves 
to demonstrate the correct behavior of the two conservation con- 
stants. For simplicity ^consider the plasma to be stationary in the 
laboratory frame (i.e., V g — 0). For circularly polarized waves. Ye 
and are constants given, in the space- independent frame, by^ ' 
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where v = e ^niax^ m e coc as defined, in Eq. (25). Substitution of Eq . 
(46) into Eqs . (43) and (45), respectively, yields 


D 2 = -n(l - n 2 ) Cl + y 2 v 2 )^ + y(l + , 

W 2 = v 2 / 2 + (1 - Ti 2 )~ i5 [(l + v 2 )" 1 +• (1 + y 2 v 2 )*Vy] 
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Note that D 2 i. = 0 and' (du e /dr 2 ) = v for circularly polarized waves. 
Eqs. (47) and (48) show clearly that the energy-momentum conservation 
parameters are functions of the wave velocity, the wave amplitude and 
the electron- ion mass ratio. In addition, from Eqs. (47) and (48) 
one obtains the following relation 


D 2 = yn(v 


72 - 


tf 2 )z 


(49) 


which shows explicitly that for a given amplitude parameter V? 2 there 
is only one value of the dynamic parameter D 2 that satisfies self— . 
consistent solutions. 


CONCLUSION 

It has been shown that the equations governing the propagation 
of relativistic electromagnetic waves in an electron-ion plasma can 


be reduced to two energy-momentum conservation relations. The two 
conservation constants, the amplitude parameter and the dynamic 
parameter, respectively , are functions of the plasma stream velocity 
the wave velocity, the electron- ion mass ratio and the wave in- 
tensity. In order to obtain self-consistent solutions the. dynamic 
parameter must be chosen appropriately when other parameters have 
been specified. 
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